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ABSTRACT 

We report the likely detection of near-infrared 2.29 /im first overtone Carbon 
Monoxide (CO) emission from the young supernova remnant Cassiopeia A (Cas 
A). The continuum-subtracted CO filter map reveals CO knots within the ejecta- 
rich reverse shock. We compare the first overtone CO emission with that found 
in the well-studied supernova, SN 1987A and find ~30 times less CO in Cas A. 
The presence of CO suggests that molecule mixing is small in the SN ejecta and 
that astrochemical processes and molecule formation may continue at least ~ 300 
years after the initial explosion. 



Introduction 



Supernovae (SNe) are suggested to be the first molecular factories in the early Universe 
( ICherchneff fc Lilly II2008I ). Molecules such as carbon monoxide (CO) are sensitive diagnos- 
tics of the temperature, density, and degree of mixing in the SN ejecta. Understanding the 
chemistry within the ejecta, and the abundance and mixing of CO and SiO molecules, is 
particularly important in studying dust formation processes since this controls how much 
carbon is available to form amorphous C grains and how much oxygen available beyond 
that bound in CO and SiO molecules to form MgSiOs and Mg2Si04. The existence of CO 
in supernova ejecta suggests that dust formation can occur, since molecules are indicative 
of relatively cool eject a (typically few thousand degrees). CO itself is an effective coolant 
( iLiu &: Dalgarnd Il995l ). so that onc e it forms it could enha nce dust formation. Isotopic 
compositions of meteoric dust grains (jClayton fc Nittlerll2004l ) show that carbon-rich dust is 
condensed within the interior of supernova ejecta. In 0-rich ejecta, radioactivity and mixing 
processes produce 'free' carbon in which to form dust. 

The quantity of dust formed within SNe ejecta has become a subje ct of debate in recent 
years (see Rho et al. 2008 and references therein; Dunne et al. 2008). iLagage et al.l (119961 ) 
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first showed evidence, in Infrared Space Observato ry (ISO) images, for dust formation in the 
nearby, young supernova remnant (SNR) Cas A. iRho et aLl (120081 ) identified three distinct 
ejecta dust populations based on their continuum shapes in Spitzer Space Telescope spectra, 
including broad features peaking at 21/im and estimated a total dust mass of between 0.02 
to 0.054 Mq. These observations provide definitive evidence of freshly-formed dust in a 
young supernova remnant and show that SNe are important dust c ontributors to the large 



quantities of dust observed at high redshifts. Theoretical models (ITodini fc Ferraral 12001 



Nozawa et al.l l2003l : iBianchi fc Schneiderl 120071 ) predict a dust mass of 0.3-1 Mq per SN. 
However, our current understanding of dust chemistry and composition is limited. The 
current models assume that very little carbon is locked up in CO and almost all of carbon 
is available for dust formation. The results from dust evolution models therefore depend 
significantly on the C/0 ratio in the gas (see Morgan & Edmunds 2003 and Dwek et al. 
1998). 

Images obtained with the Spitzer Infra red Array Camera (IRAC) each present a different 
view of the nucleosynthesis within Cas A (lEnnis et al.ll2006l ). A surprising result from the 
IRAC observations is the similarity between the morphology of the IRAC Band 2 image 
(centered at 4.5 /im) and ejecta-dominated images of the remnant. The origin of this emission 
is not known, bu t based on the lack of bright Ha emission, the contribution from Bra emission 
is likely small. lEnnis et al.l (120061 ) therefore suspected that CO emission was responsible 
for the emission at this wavelength, but could not be confirmed since the Spitzer Infrared 
Spectrograph (IRS) has a blue cutoff of 5.2 /im. 

As a tracer of gas properties and dust formation, CO is potentially an important di- 
agnostic for any supernova. The first detectio n of CO from supernova ejecta was from 



the well-studied. Type II supernova, SN1987A (ICatchpole et al.l Il989l : iMeikle et al.l Il989 



Spyromilio et al.lll988l ). providing an unique insight into the chemical dynamics of the ex- 
plosion and ejecta. The variation of the observ ed CO emission with time after the inital 
explosion was modeled by lSpyromilio et al.l (119881 ) . who suggested that CO formed out in the 
SN ejecta at a temperature of 1800 K and expansion velocit y 1200 km s~^ at 337 days af t er the 



explosion. They determined a total mass Mco ~ 10 ^Mq. ISpyromilio fc LeibundgutI (119961 ) 
detected ~ lO'^M© of CO forming at temp erature 4000 K in SN 1995ad. CO has further 



been detected from the Type Iln S N 1998S (ICerardv et al.ll2000l) . SN 2002dh (IPozzo et al. 



2006h . and the Type Ic SN 2000ew JCerardv et al.ll2002h . The number of CO detections m 
SNe in the literature is still relatively small, and all are restricted to within a few years after 
the initial explosion. In this Letter, we report the likely detection of the CO overtone from 
Cas A, an approximately 300-yr old supernova remnant. This is the first time CO emission 
has been detected in a Galactic supernova remnant and suggests that molecule formation is 
not only a common occurrence in core-collapse ejecta but molecules also exist (or continue 
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forming) long after the initial explosion event. 



2. Observations 

We observed Cas A on 2006, September 3 and 4 with the Wide field InfraRed Camera 
(WIRC, Wilson et al. 2003) on the Hale 200 inch (5 m) telescope at Mount Palomar. The 
weather was very clear with the seeing typically 0.8". The camera field of view is 8.7'. We 
observed a set of dithered positions (designed to uniformly cover the supernova remnant), 
with interleaved, dithered observations of an 'off' position 30' away. A sky reference image 
was generated from the 'off' images, and was subtracted from the 'on' images after scaling 
to match the median sky brightness. This observation and analysis procedure preserved the 
diffuse emission of the SNR. We took narrow-band images using filters centered on CO (2-0) 
at 2.294 /xm, K-band and H-band continuum at 2.270 and 1.570 //m, and P/3 at 1.182 fim 
(see Figured]). The narrow-band filter widths are typically 0.025 ^m. The exposure time of 
CO, K-cont and P/9 is a combined 90 sec from three back-to-back images of 30 sec exposure 
(to reduce the number of saturated stars), with a total on-source integration time of 8100 sec 
for CO, 7290 sec for K-cont and 810 sec for P(3 image. The exposure time of H-cont is 90 sec, 
with a total on-source integration time of 1440 sec. Respective sensitivities of the narrow- 
filter images are 1.70x10"^ erg s~^ cm~^ sr~^ for CO, 1.04x10"^ for K-cont, 1.54x10"^ for 
H-cont, and 1.25 xlO~^ for P/3. T he detailed data reduction methods are similar to those 



described in lKeohane et al.l (120071 ): additionally, focal plane distortion was derived from the 
images and corrected. We directly subtracted the K-continuum image from the CO filter 
i mage in order to subtract stars and contamination from synchrotron emission from Cas A 



( iRho et al.ll2003l ). Due to excellent seeing (0.8") and low airmass towards Cas A, the stars 
were well subtracted. We further masked saturated stars and applied adaptive smoothing 
and wavelet filtering to the images. 



3. Results 

Figure [T] shows a color composite image of Cas A, combining the P/? line (blue), 2.27 
/im continuum (green), and CO (2-0) 2.29 //m line (red). P/3 traces recombining atomic 
hydrogen. The 2.27 //m continuum consists mostly of non-thermal synchrotron radiation 
and delineates the shell structure (see Rho et al. 2003). Finally, the CO (2-0) overtone 
emission traces the dense knots of gas along the shell and interior interface. The molecular 
emission is evident as red in Figure [T] and more easily seen in the continuum-subtracted CO 
image. Figure [H revealing the line emission towards the northern parts of bright rings and 
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an eastern filament of tlie Minkowski knots in tlie soutlieast. Tlie CO emission is liiglily 
clumpy and distributed in a collection of knots with an upper limit on t he size of ^ 0.01 3 
pc (limited by the spatial resolution of 0.8") at a distance 3.41qj kpc (IReed et al.l Il995l ). 
The flux densities for the dense gas knots are listed in Table 1. We selected the brightest 
knots of 18 regions in the CO, K-continuum or P/3 images with angular sizes 3"-15" and 
estimated the flux densities after background subtraction. The regions include knots which 
are bright in the continuum-subtracted CO image and the flux densities of individual knots 
range from 0.1 to 0.3 mJy. The flux densities of two representative CO knots and one other 
knot in the NW are given in Table 1 along with the positions. The total flux density of 
the CO-detected area after K-continuum subtraction (see Fig. 2a and 2c), is 1.27±0.17 
mJy, where the estimated error includes calibration (7%), background variation (5%) and 
statistical (10%) errors. Some faint, diffuse CO emission from the Cas A may have been 
removed from our image because of the continuum subtraction; the flux reported here is 
conse rvative, including only th e bright knots. After ext inction correction using A„ of 4.5 
mag (IFesen &: Gundersonlll996l ) and the extinction law of iRieke &: LebofskyI (119851 ). the flux 
density is 2.03±0.27 mJy, equivalent to a CO overtone luminosity of 0.01±0.0014 d3.4fcpc 
Lq. The ratio of H-continuum to K-continuum is globally consistent with the origin of 
synchrotron radiation with a synchrotron index of 0.65 - 0.71 (Rho et a l. 2003); moreover, 
the H-continuum and K-continuum images are similar to the radio image (jWright et al.lll999 
Rho et al.ll2003l ). indicating they are dominated by synchrotron emission. 



We compiled spectral energy distributions (SEDs) of near-infrared H and K bands and 
Spitzer IRAC Bands 1 to 3 (3.6 - 5.8 i^im) bands for 18 knots; Figure [3] shows three repre- 
sentative SEDs after correcting for the extinction. The SEDs show a sudden jump at 2.294 
/im toward the CO detection positions, with no such jump toward the NW ejecta position. 
The SEDs of the CO knots also show a larger IRAC Band 2 (4.5 /im) excess than that of 
the NW ejecta knot. 

A few knots towards P(3 emitting regions (blue in Figure 1) show K-continuum brighter 
than the CO-band emission. These knots look like elliptical spheroids, indicative of Quasi 
Stationary Flocculi (Fesen & Gunderson 1996, Willingale et al. 2002). We further note that 
H-continuum emission towards hydrogen knots is slightly higher than that expected from 
synchrotron emission, suggesting that despite the narrow bandwidth of H-continuum filter, 
the image may be picking up the weak hydrogen line Br A at 1.57 /ir n. The CO filter d o es not 
include any hydrogen or helium lines. The near-infrared spectra of iGerardy fc Feseru (120011 ) 
did not detect any hydrogen and helium lines within the H-cont, K-cont, CO filters, but our 
images have a higher sensitivity and therefore may be detecting such lines. 



The P/3 emission is largely from shocked circumstellar gas, and similar to the optical 
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Ha emission (IFesen fc Gundersoru Il996l : iFesenI 1200 ll ) . The emission is mainly clustered in 
the northern-most inner shell and outside the southwestern shell, as shown in Figure [H 
The former contains both circumstellar gas and fast moving knots. The hydrogen lines of 
circumstellar medium knots als o show moderate high veloc ity components, which are defined 
as "mixed emitting knots" by iFesen fc GundersonI (119961 ). who discussed their origin. In 
particular, rich infrared ejecta e mission is detected f r om th e northern most inner shell (see 



Figure 2 of Rho et al. 2008). iFesen fc GundersonI (119961 ) suggested that mixed emitting 



knots may be "splatter" fragments produced as high speed underlying core condensations 
broke through and partially mixed with the progenitor star's surface layers. 



4. Discussion 



We detected significant near-infrared CO overtone emission from small portions of the 
Gas A remnant towards the northern part of the bright ring and in the eastern region of 
Minkowski knots in the southeast. It is possible that we are missing fainter GO emission 
elsewhere in the SNR due to our detection limit (of 1.98x10^^ erg s~^ cm~^ sr~^) in the 
continuum-subtracted GO ima ge. The locations where GO emission is detected coincide with 



the eiecta at the reverse shock (IHwang et al.ll2004J : iFesen fc Gundersonlll996l : iRho et al.ll2008 
Yang et al.ll2008l ). It indicates that the GO gas was formed in the ejecta. 



The Spitzer IRAG Band 2 image shows a morphology similar to ejecta-dominated im- 
ages, but the origin of this emission was a mystery. After scaling from the overtone to 
fundamental band, our observations suggest that the GO fundamental contributes some 
(>l/6) but possibly not all of the emission in IRAG Band 2 (Figure [3]). Other contributors 
to this band include multiple [Fe II], hydrogen recombination and helium lines. The ratios 
of IRAG 2 and GO-filter emission from 18 knots across Gas A suggest that GO contributes 
10-30% of the emission in the GO emitting regions when compared with non-GO emitting 
ejecta knots. Gonstraining the percentage of GO contribution to the IRAG band 2 image 
is difficult since the ratio of the fundamental to the first overtone varies depending on the 
ejecta temperature and density. We cannot therefore strongly constrain the GO contribution 
to the IRAG Band 2. We analyzed the IRS spectroscopy cube to find any signature of the 
fundamental GO (4.8 - 6 /im), but the signal-to-noise was too poor. The Spitzer IRS data 
are not sensitive enough to detect this emission (errors of the IRS spectra were 1-2 MJy sr~^, 
an order of magnitude higher than the typical fiux densities of near-infrared GO emission) 
although there might be some hint for the northern position (at the position of the strongest 
21 /im dust spectrum). We also looked for SiO molecules (at 8-10 /im), but could not 
distinguish them from the silicate emission feature at 9.8 /im (silicate stretching mode). 
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We estimated the mass of CO in Cas A by comparing the CO flux density in Cas A 
with that found in the supernova SN 1987A, since without spectra it is not possible to derive 
a density, temperature, and velocity of CO gas. Figure B | comp ares the CO emission from 
Cas A and SN 1987A jMeikle et all llQsd : IWooden et all Il993h . After accounting for the 
distance differences, the CO (2-0) first overtone flux density of Cas A is a factor of ~ 28 
smaller than that of SN1987A on day 255. Assuming the same gaseous conditions as seen 
in the CO emitting regions in SN1987A (a temperature of 1800-2800 K and a velocity of 
2000 km s~^), we estimated the CO mass in Cas A to be Mco ~ 1.7 x 10~^Mq for Local 
Thermodynamic Equilibriura (LTE) (jSpyromilio et al.lll988l ) and 8.2 x 10^^ Mq for non-LTE 
( iLiu. Dalgarno fc Leppl[l99i ). Note that the CO mass estimate is highly uncertain because 
we assume the gas properties described above. 

The CO detections in SNe have been within a few years after the i nitial explosion. The 
mairi destruction process of CO is from impacts with energetic electrons (ILiu. Dalgarno fc Lepp 
19921 ) and with ionized gas atoms. The latter process quickly destroys CO, hence the de- 
tection of the CO fundamental band provides constraints on the degree of mixing in the SN 
ejecta. Given that Cas A is believed to be the remnant of a Type lib explosion (Krause et 
al. 2008), this suggests that either the CO layer is not macroscopically mixed with ionised 
helium or that the helium is not ionised. Typical elect ron density of C- and 0-rich region s 
during formation of CO is shown to be 5x10*^ cm~^ ( IGearhart. Wheeler fc Swarta Il999l ). 
The highly clumpy structures can be explained by such high density CO gas within the 
warm or hot ejecta, which naturally fortifies t he evidence of fresh dust formation observed 
in Cas A jLagage et al.lll996l : iRho et al.l[2008[ ). 



In dust evolution models, understanding CO molecules in SN ejecta is a key physi- 
cal question. The formation efficiency of CO molecules critically affects the dust forma- 
tion efficiency. In chemical equilibrium, the abundance of CO is balanced between radia- 
tive processes and radioactive decay fo ^^Co, because the main de struction process is by 
the energetic electrons produced by the radioactive decay of ^^Co (ILiu. Dalgarno fc Lepp 



19921: IClavton. Deneault fc Mever Il200lh . Both Nozawa et al. (2003) and iTodini fc Ferrara 
2001) assumed the mass fraction of C atoms that are not locked in CO is 99-100%, and 



Travagho et al.l fll999l ) assumed that mixing occurs at the molecular level, prior to dust con- 
densation. In contrast, C l ayton & Nittler (2004), Deneault, Clayton & Heger (2003), and 
Fryxell. Arnett fc Muellerl (jl99ll ) suggest that gases are unlikely to be mixed at the molec- 
ular level within a few years. Our CO detection suggests that molecular mixing is small (at 
least smaller than previously thought) in the SN ejecta and during the development of the 
reverse shock. CO molecules in Cas A suggests a higher CO abundance in the gas phase and 
lower efficiency in dust formation of carbon dust. Remaining molecules also indicate that 
dust formation was inefficient (leaving molecular leftovers) or occurred inhomogeneously (e.g. 
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only in dense knots). Deep infrared spectroscopic observations will allow us to measure the 
velocity, temperature and density of the CO molecules in Cas A and to study the on-going 
chemistry of molecules in SN ejecta. 

5. Conclusions 

We detected carbon monoxide emission from the young supernova remnant Cas A. This 
may be the first CO detection from an older SN (age ~300 yr). CO emission has also been 
detected in 8 other core-collapse remnants, but in all cases it has been limited to within 
3 years of the explosion. Theoretical models suggest that dust formation in SNe occurs 
between 200-800 days after the SN event yet the observations presented here suggest that 
molecular and dust formation processes could continue long after. The detection of CO in 
the Cas A remnant demonstrates that astrochemical processes are still on-going more than 
300 years after the initial explosion. 
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comments. H. Gomez would like to thank Las Cumbres Observatory for their support. 
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Based on observations obtained at the Hale Telescope, Palomar Observatory as part of a 
continuing collaboration between the California Institute of Technology, NASA/JPL, and 
Cornell University. 



Table 1: Summary of Observed Brightnesses 



filter 


flux density" 


Flux Density (m 


Jy) 




range 


CO knotl" 


CO ^012" 


NW knot'' 


2.27 iim (K-cont) 


0.1-1.0 


0.610±0.009 


0.290±0.006 


1.780±0.016 


2.29 fim (CO) 


0.2-1.2 


0.850±0.013 


0.490±0.010 


1.880±0.020 


3.6 Aim (IRAC 1) 


0.4-2.3 


0.71±0.20 


0.38±0.16 


2.23±0.35 


4.5 Atm (IRAC 2) 


1.2 -15.0 


11.03±0.64 


5.10±0.45 


9.24±0.60 


5.8 iim (IRAC 3) 


0.6-4.2 


2.20±0.42 


1.04±0.34 


3.05±0.67 



"The range of flux densities (in mjy) towards CO detected positions. ''Positional coordinates in RA and Dec (J2000): CO knotl 
(green in Figure[3ll at 350.850°, 58.840°, knot2 (red) at 350.854°, 58.834°, and the northwest ejecta knot (blue) at 350.825° and 
58.836°. The flux densities were measured by averaging over an angular size of 5.3"x3.1", 38"xl0", and 6"xl2" for knotl, 
knot2, and the NW knot, respectively. 



- 8 - 




Fig. 1. — Mosaicked color image of Cas A with WIRC, with CO (red), K-continuum (green) 
and P/3 (blue). The range of surface brightnesses measured from the diffuse structures are 
0.31 - 2.79x10"^ for CO, 0.23 - 1.4x10"^ for K-continuum, and 1.2 - 29.0x10"^ erg s"^ cm'^ 
sr~^ for Pj3. CO-excess regions (in red) are marked as boxes (the bottom box marks the 
location of the Minkowski knot). The image is centered at R.A. 350.858° and Dec. 58.814° 
(J2000), and covers an 7'x7' arcmin field of view. 
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Fig. 2. — Zoomed-in continuum-subtracted-CO (a and c) and IRAC2 (b and d) color images 
towards CO-detected regions. These zoomed-in images cover the regions marked as boxes in 
Figure [D the pannels (a) and (b) show the northern bright ring [centered at R.A. 350.857° 
and Dec. 58.836° (J2000), with area 2. 5' x 1.2'], and pannels (c) and (d) are towards the 
Minkowski knots [centered at R.A. 350.888° and Dec. 58.802° (J2000), with area 1.2'x0.6' ], 
respectively. The near-infrared CO images (a and c) include the CO (2-0) first overtone, and 
the Spitzer IRAC 2 images (b and d) include CO (1-0) fundamental. The locations of the 
three knots in Table 1 are marked as green rectangles (upper: CO knotl, lower: CO knot2, 
and right: NW knot). 
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Fig. 3. — The spectral energy distribution of two representative CO emitting knots: knotl 
(green) and knot2 (red) (positions are given in Table 1) which show excess emission at 
2.294 fim. The contribution from the synchrotron continuum is shown by the dotted line. 
Comparison with the SED (scaled to the CO knotl continuum) of the NW ejecta knot (blue 
dashed curve) shows that the knots with near-infrared CO-filter excess emission may also 
contain the CO fundamental emission in the IRAC band 2. The error bars of near-infrared 
emission include both statistical and calibration errors. The SN 1987A spectrum is shown 
for comparison (orange curve) at 255-260 days after the explosion event (courtesy of Peter 
Meikle, Meikle et al. 1989; Wooden et al. 1993). 
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